Implementation of a Quantum Annealing Algorithm Using a Superconducting Circuit 
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A circuit consisting of a network of coupled compound Josephson junction rf-SQUID flux qubits 
has been used to implement an adiabatic quantum optimization algorithm. It is shown that detailed 
knowledge of the magnitude of the persistent current as a function of annealing parameters is key 
to implementation of the algorithm on this particular type of hardware. Experimental results 
contrasting two annealing protocols, one with and one without active compensation for the growth 
of the qubit persistent current during annealing, are presented in order to illustrate this point. 

PACS numbers: 85.25.Dq, 03.67.Lx 



The successful implementation of any solid state quan- 
tum information processor will ultimately depend upon 
having defeated several critical challenges. While noise 
[1] and fabrication variability [2] arguably gate progress 
at the moment, it must be recognized that practical issues 
related to scalability, architecture and algorithms also de- 
serve attention. It may be within experimental grasp to 
address some of these issues in the context of adiabatic 
quantum computation [3, 4] using state of the art designs 
and fabrication methods. In this article, we address sev- 
eral practical details concerning the implementation of an 
adiabatic quantum optimization algorithm using a net- 
work of coupled rf-SQUID flux qubits [5]. It is shown 
that detailed knowledge of qubit properties as a function 
of annealing parameters is key to successful implemen- 
tation. Experimental results from a chain of six coupled 
qubits subjected to two different annealing protocols are 
presented in order to illustrate this latter point. 

The work presented herein focuses on hardware de- 
signed to enable a particular adiabatic quantum opti- 
mization algorithm [4] for computing the vector s that 
minimizes the objective function 
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where si = ±1 and N is the length of s. Here, hi and Kij 
are dimensionless real numbers that arise from a partic- 
ular choice of problem instance. This type of problem is 
of interest as it is known to be NP-hard [6]. Equation 
(1) can be recast as the potential energy of a system of 
coupled spin-1/2 particles via the substitution s t — > a* , 

z(x) 

eigenstates of a z be denoted by ||) and ||). The vector 
s that minimizes Eq. (1) is then encoded in the ground- 
state |<7f) = |si . . . sn) of an Ising spin glass. The algo- 
rithm for finding Iqf) relies upon exploiting the trans- 



glass. Let the Hamiltonian of such a system be 

H(t) = A(t)Hi + B{t)H F (2) 
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where 0<t<t f , A(0)/B(0) > 1 and A(t f )/B(t f ) < 1. 
We refer to A(t) and B(t), which have units of energy, as 
envelope functions. If A(0) is much larger than all other 
relevant energy scales, including temperature, then the 
system will begin the evolution at t = in the ground- 
state of Hi, \gi) = nili (II) + IT))/V2, with a proba- 
bility of 1. If the subsequent evolution is adiabatic, then 
the system will be found in \gp) at t = tt. Note that 
the ^-dependence of this algorithm is entirely contained 
in A{t) and B(t), which are explicitly not functions of N 
or problem instance. 

The most convenient forms for the envelope functions 
depend upon the details of the hardware. Consider a net- 
work of coupled compound Josephson junction (CJJ) rf- 
SQUID qubits [2, 7] whose persistent currents 
and tunneling energies A($^j) are identical as a func- 
If is a function of t, then the 



tion of CJJ bias $* {{ . 



processor Hamiltonian can be expressed as 
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where €i(t) = 2 3>f (t) represents the energy bias of 

a flux qubit (-$ < *3j(f) < -$o/2, $o = h/2e), $f (t) 

i i ° 

is an externally controlled flux bias and J^- = My \Ig(t) | 



where a\' {x) is the z(x) Pauli matrix for spin i. Let the represents the pairwise coupling mediated by a mutual 

inductance My. Comparison of Eqns. (2) and (3) read- 
ily yields one envelope function: A(t) = A(t)/2. On the 
other hand, there is no unique definition of B(t). One ap- 
proach is to scale €i(t) and Jy(i) by a convenient factor: 
verse (jJx') degrees of freedom of a quantum Ising spin let B(t) = Jafm(^) = Mafm Uo(^)| 2 , where Mafm is 
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FIG. 1: (Color online) Optical image of circuit with qubits 
qi . . . Q6 labeled. Inset shows an electron micrograph of 
with bias lines for controlling qubit CJJ bias (5?cjj)> qubit flux 
bias ( &q), readout flux bias (<&f ) and readout current bias 
(i ro ) as indicated. Transformers between qubit and couplers 
denoted as M co - q . Transformer between qubit and readout 
denoted as M ro -„. 



the strongest antiferromagnetic (AFM) coupling needed 
to embed a particular problem instance. Doing so implies 
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which yields a prescription for mapping Eq. (1) onto 
the hardware: h, = ^(^/Mafm \lq(t)\ an d ATjj = 
Mij /Mafm- Note that ify has no i-dependence. To 
avoid i-depcndcncc in hi one must use a i-dependent 
qubit flux bias: 



^(0 = ^m afm |4 p W| 



(5) 



We denote annealing processes that use $f (t) as defined 
by Eq. (5) as controlled annealing. Processes in which $f 
are held static will be termed fixed annealing. Formally, 
the relative contributions of ai - versus Oz^o* -terms 
to Eq. (3) will vary during fixed annealing. It will be 
demonstrated that fixed annealing is not a viable means 
of implementing an optimization algorithm and that con- 
trolled annealing remedies the problem cited above. 

We have performed a series of annealing experiments 
on a chain of 6 CJJ rf-SQUID qubits with intervening 
AFM couplers. These particular qubits, as depicted in 
Fig. 1 and labeled as qi . . .qe, were part of a larger cir- 
cuit and were selected on account of their low Joseph- 
son junction asymmetry [2]. The chip was fabricated 
on an oxidized Si wafer with Nb/Al/A^Oa/Nb trilayer 
junctions and three Nb wiring layers separated by sput- 
tered Si02- It was mounted to the mixing chamber of 
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FIG. 2: (Color online) CJJ bias dependence of |/|| and en- 
velope functions A = A/2 and B = Japm = Mafm \l$ \ 2 ■ 
Direction of CJJ bias ramp during annealing as indicated. 



a dilution refrigerator and cooled to T = 35 mK in the 
presence of a very low (< 9nT) background magnetic 
field inside a PbSn coated shield. Each qubit was con- 
nected to three others via in-situ tunable rf-SQUID cou- 
plers, which we treat as classical effective mutual induc- 
tances [8]. Each qubit was also inductively coupled to a 
hysteretic dc-SQUID for readout [9] The chain of qubits 
studied herein was isolated from the rest of the chip by 
tuning unused couplers to provide zero coupling and bias- 
ing unused qubits with = — < &o/2 to minimize their 
persistent currents. The CJJ bias dependences of 
and A of these qubits were synchronized by applying the 
methods described in Ref. [2]. Further details regard- 
ing qubit calibration and measurements of (/^(^Jjj) | and 
A($^jj) can be found therein. The expected CJJ bias de- 
pendence of the envelope functions, as determined from 
the mean device parameters reported in Rcf . [2] , has been 
plotted in Fig. 2. Here, the synchronization CJJ bias 



has been defined such that A($^)/2h = 10 MHz. 
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FIG. 3: (Color online) AFM coupled domain experiments. 
Qubits represented as spin- 1/2 particles (circles with arrows) 
and AFM couplers as solid lines, n denotes number of qubits 
per AFM domain (dashed boxes). Flux biases 3>f as noted. 
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FIG. 4: (Color online) Waveforms for a) fixed and b) con- 
trolled annealing. All active qubits utilize the same wave- 
form. Only qubits q$ and qi carried i-dependent $f . 



To study controlled versus fixed annealing, we have 
measured three different configurations of the chain, as 
shown schematically in Fig. 3. The diagram depicts a se- 
quence of experiments in which qubits and were sub- 
jected to variable flux biases while all other qubit flux bi- 
ases $f = 0. Qubits that were not used in a particular ex- 
periment had their CJJ biases held at = — $o/2 so as 
to decouple them from the active qubits. In each succes- 
sive experiment, one more AFM coupled qubit was acti- 
vated on both ends of the chain. In the limit A <C Jafm, 
one can write analytical solutions for the eigenstates of 
Eq. (3) for each of the experiments depicted in Fig. 3. 
The results show that the four lowest energy levels can 
be ascribed to a system of two AFM coupled effective 
qubits with the same |/^| as a single qubit, but a renor- 
malized tunneling energy A ff ss A"/(2Jafm)' 1_1 , where 
?i is the number of qubits in an AFM domain: 



H(n, t) » - £ 1 [ £i (i)cr« + A cff (n, t)a% 

i=3 
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Thus, the experiments depicted in Fig. 3 were isomorphic 
2-qubit optimization problems in which the tunneling en- 
ergy A e g was a strong function of n. An optimization 
algorithm applied to these systems ought to yield results 
(states of (?3 and 94) that are independent of n. 

The fixed and controlled annealing waveform patterns 
are depicted in Fig. 4a and 4b, respectively. In both 
cases, the CJJ biases of all active qubits were simul- 
taneously ramped from — $o/2 to —$0- The ramps 
were digitally low pass filtered using f c = 1.3 MHz so 
as to avoid uncontrolled delays between qubits due to 
the limited bandpass of the wiring [2] . The slope of the 
ramps in the vicinity of $^j(to) = was ^cjj/^ = 
3.2 x 10 -2 3>o/ms. For fixed annealing, $3 and $| were 
set to static values during the CJJ ramps. For controlled 
annealing, $f(t) and 3>f(£) carried scaled time depen- 
dent waveforms as dictated by Eq. (5). Here, 
was obtained by sampling the smooth function shown 
in Fig. 2. Choosing the target values of = £3/2Jafm 
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FIG. 5: (Color online) Most probable state map | (5 1 53 94)! 2 
for the n = 1 experiment in the limit A <C Jafm- Dashed 
lines indicate phase boundaries in the limit A — » 0. 



and = £4/2 Jafm and knowing A/afm = 1.35±0.02 pH 
from independent measurements [2] allowed for construc- 
tion of <f>! (t) and $!(£)• In both fixed and controlled 
annealing, these fluxes were set to zero during qubit ini- 
tialization and at the end of annealing prior to readout. 

To develop a sense of how A e g(n) impacts the dy- 
namics, consider the n = 1 experiment as a function of 
£3/2 Jafm and £4/2 Jafm in the limit A <C Jafm- The re- 
sults of a calculation of the probability | (<? | (73*74) | 2 , where 
\g) is the groundstate of Eq. (6), are shown in Fig. 5. 
Here, the AFM states \\f\) an d IT!) occupy the majority 
of the plane and intersect at a diagonal phase boundary. 
Horizontal and vertical phase boundaries are located at 
±1 where the AFM regions intersect the ferromagnetic 
(FM) states |J.j) and | TT> - The width of the region about 
the phase boundaries where one may observe quantum 
tunneling between states are as noted in the diagram. 

Let the system depicted in Fig. 5 be subjected to fixed 
annealing. In this case, ei(t)/2J\py[(t) cx 1/ 
Since is a monotonically increasing function of t, 
then (<!>! , $f) corresponds to a point that moves radially 
inward in the (£3/2 Jafm, £4/2 JAFM)-plane as annealing 
progresses. If this point traverses a phase boundary at a 
time t p , then it will be forced through a phase transition. 
The probability of faithfully tracking the groundstate will 
be a function of A(i p ) and of the time spent in the vicin- 
ity of the phase boundary, which will be proportional to 
d&cjj/dt\t=t p - In this regard, fixed annealing could be a 
generalized Landau-Zener experiment [10]. Now consider 
replacing the single qubits (73 and (74 with AFM domains 
of size n > 1: €i{t)/2 Jafm(£) will remain unchanged, 
but A(i) — * Aoff(i). Since the tunneling energy will be 
suppressed, the probability of achieving the groundstate 
will decrease with n when all other experimental param- 
eters are equal. Thus, the locations of apparent phase 
boundaries in the (<&f , C^-plane will depend upon n. 

In contrast to the above scenario, a system subjected to 
controlled annealing never traverses a phase boundary in 
the limit A <C Jafm- This is assured because ei/2J^FM 
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FIG. 6: (Color online) Maps of the most probable final spin 
configuration as a function of fixed annealing flux biases (a, 
b and c) and controlled annealing dimensionless parameters 
(d, e and f). States are labeled as \q3q4) (a and d), \q2q3q4q5) 
(b and e) and \q1q2q3q4q5q6) (c and f). 

will be independent of t by construction. As such, the 
locations of phase boundaries in the controlled annealing 
(e3/2Jafm, £4/2 Jafm) -plane will be independent of n. 

Experimental maps of the most probable final spin con- 
figuration for fixed and controlled annealing were gener- 
ated by sampling on a grid of points in the ($3,^4)- 
and (£3/2 Jafm, £4/2 JAFM)-plane, respectively. Since the 
state of each qubit could be read by its own dedicated 
dc-SQUID magnetometer, it was possible to unambigu- 
ously identify the final spin configuration of any given 
AFM domain at the end of every annealing cycle. Run- 
ning 64 cycles per point provided histograms of the final 
spin configuration from which one could readily identify 
the most probable state. We note that there were no 
measurements that ever indicated single qubit flips had 
occurred in any of the AFM domains. 

Results of the fixed annealing experiments are shown 
in Figs. 6a^c. Here, we have colored each point in 
the (<I>|, $|)-planc to represent the spin configuration 



I - • - ?3 94 - ■ • ) that was observed with the highest proba- 
bility. The data indicate that the apparent AFM interac- 
tion between q 3 and (74 becomes weaker with increasing n, 
as evidenced by the progressive encroachment of the FM 
states. The difference in flux between the pairs horizontal 
and vertical boundaries can be interpreted as an effective 
AFM coupled flux 2$* = 4M AFM where ^ 

represents a particular CJJ bias that depends upon n. 
Decreasing <£* indicates that the dynamics of the AFM 
chain effectively freeze out at larger where is 
smaller (see Fig. 2). A detailed study of the dependence 
of 4>* upon n and d&^/dt may provide information re- 
garding the mechanism by which this circuit achieves its 
final configuration. This matter will be the topic of a 
future publication. Nonetheless, these data clearly show 
that fixed annealing is not a viable means of operating 
this hardware as a processor as the solution to a given 
isomorphic 2-qubit problem depended upon n. 

The controlled annealing results presented in 
Figs. 6d— >f show no dependence upon n. To within 
experimental error, the phase boundaries in the 
(£3/2 Jafm, £4/2 JAFivi)-plane agree with the results of 
Eq. 6 in the limit g — » 0. Additional measurements also 
revealed no dependence upon d$^/dt over 3 orders of 
magnitude (3.2 x 1CT 4 to 0.32$ /ms). Therefore, it 
has been demonstrated that the controlled annealing 
protocol returns the correct solutions to optimization 
problems posed as a set of hi = 6^/2 Jafm and = lj 
per Eq. (1). Consequently, this circuit, when used in 
conjunction with controlled annealing, can be viewed as 
a prototype quantum Ising spin glass computer. 

Conclusions: A method for annealing coupled CJJ rf- 
SQUID flux qubits that accounts for realistic device be- 
havior and embodies the physics of a particular quan- 
tum adiabatic optimization algorithm has been experi- 
mentally demonstrated. This work represents a critical 
step toward the development of practical adiabatic quan- 
tum information processors. 
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